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MEUTROh1  DIFFUSION  AT  LARGE  DIST/f’CES 

fro;  the  source. 


by 

K . T . 3 p I n n e y 


Summary. 


k theory  of  neutron  diffusion  by  11.  A.  Be  the  is  generalised 
to  include  the  case  of  capture  and  a nixed  medium. 
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1.  Introduction 


Bethe^  has  provided  a theory  which  deals  with  neutron  densities 
at  large  distances  from  the  source  and  covers  the  "gap”  between  age 
theory  and  the  asymptotic  theory  of  Trick.  This  theory  has  obvious 
applications  for  calculating  the  efficiency  of  neutron  shields,  except 
that,  as  presented  in  the  above  paper,  it  considers  only  a single  type 
of  nucleus  in  the  scattering  medium  and  takes  no  account  of  capture. 
However  it  may  be  easily  generalized  to  include  the  case  of  a mixed 
medium  and  constant  capture  and  it  is  the  purpose  of  this  paper  to  do 
so. 


It  is  assumed  that  the  reader  is  familiar  with  Bethe’s  paper  and 
the  argument  is  not  elaborated  in  detail  wherever  it  follows  closely 
that  of  Bethe.  We  use  the  same  notation  except  where  differences  are 
noted  in  the  text. 

All  scattering  and  capture  cross  sections  arc  taken  independent 
of  neutron  energy  and  we  consider  only  elastic  collisions. 

As  in  Bethe’s  paper  we  shall  limit  ourselves  to  an  isotropic  mono- 
• energetic  plane  source  in  an  infinite  homogeneous  medium. 
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2.  The  Boltzman  Equation 

Suppose  the  medium  contains  nuclei  of  m different  types  with  masses 
Mj,  j = 1 ...  m and  Mj  » 1 for  all  j. 

Let  the  scattering  cross  sections  of  the  nuclei  be  er^  j ■ 1. ..  m 

and  the  capture  cross  sections  crcj  , j ■ 1 ...m 

Let  the  nuclei  be  distributed  uniformly  through  the  medium  with  Nj 
as  the  number  of  nuclei  of  type  j per  cm*. 

Then  the  total  mean  free  path,  l , is  defined  by 


1 

l 


m 

T. 

j“l 


Nj  Kj  * 


^cj) 


In  order  to  derive  the  modified  form  of  Bethe’s  eqn(l)  we  start 
from  eq”(63)  in  the  review  by  Marshak  ^ *'  on  which  he  bases  his  work. 
That  eq11  was  written  down  for  a non-capturing  medium  but  it  is  obvious 
that  it  aiso  holds  for  a capturing  medium  if  the  function  f is 
defined  properly. 

The  equation  is 
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l u — ♦ v(z,u,u) 
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where  f(uo,  u-u')  is  the  relative  probability  that  a collision  will 
change  the  logarithm  of  the  inverse  neutron  energy  from  u*  to  u with 
deflection  through  an  angle  6 given  by  uo  = cos  0.  u is  taken  to 
be  zero  at  the  source  and  v is  the  neutron  collision  density. 

In  our  case  we  must  have 


f(ooi  u— u 1 ) c 5!  Cj  fj  (no.  U-u')  ....  (2) 


where  cj  = Z Nj  o-gj, 


this  being  the  relative  probability  that  a neutron  will  be  scattered 
by  a nucleus  of  type  j in  a collision. 

For  elastic  collisions  the  functions  fj  are  given  by 
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for  u1  i u - qj 


0 for  u'  < u - qj 


where  qj 
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We  may  write  (1)  as  follows,  taking  the  total  mean  free  path  as 
our  unit  of  length 


(jw  mu  a(z)  S (u) 

u -T  + 'K(z.u.u)  = J dfl'  l J du'  v(zu'u')c;  f i ( uo,  u—u ' )* 

dz  j-1  u-qj  J J 4« 


....  (4) 


The  capture  of  neutrons  is  expressed  hy  the  fact  that  the  sum  of 
the  coefficients  Cj  is  less  than  unity. 


3.  Fourier  - Laplace  transform. 


The  convolution  theorem  for  Laplace  transforms  states 
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provided  that  u0,  as  defined  by  no  ■ e ^ — — - — 

c <5 

satisfies  0 # u0  $ qj  and  this  is  true  from  the  definition  of  qj 
maximum  possible  logarithmic  energy  less  in  one  collision. 

Now  apply  a Fourier  Laplace  transformation  to  (4)  and  we  get 
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where 


r * f00  izy  “ -nu 

<P(ynn)  = J e dz  j du  e ?(zu|i), 


Then  writing  the  last  factor  of  (6)  as  an  exponential  we  have 
(M-;  + 1' 


gj(non) 
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where  we  have 


neglected  terms  of  order  —3  and  above. 
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4.  Solution  by  expansion  in  spherical  harmonics. 

Now  expand  in  spherical  harmonics  so  that 

21  Sjj(h)  PjUo)  •••  (8) 

where 


gjitn)  - (21  + 1)  / gj(uo,n)  P^uo)  an 


etc. 

Now  expand  <P(y,nu)  in  spherical  harmonics  so  that 

<p(ynn)  - l (2 ( +i)  Pz(u)  <p»  (yti)  ...  (18) 

l 

where 

'Pt(yh)  * ^ J <P(ynu)  PjU)  dn  ...  (13) 

Then  (6)  becomes,  on  substitution  from  (8)  and  (18) 
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where  gj(n)  - l cj  gj  j (n ) 
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...  (15) 


We  now  look  for  the  value  y ■ -iv  for  which  the  homogeneous 
equation  corresponding  to  (14)  can  be  satisfied,  i.e.  we  consider 


N Pj(ii) 

<p(v,n,n)  - 2 <Pj(vn)gz(n)  

1*0  1 —v\i 


(16) 


where  we  consider  only  a finite  number  of  terms  of  the  infinite  series 
since  it  is  rapidly  convergent. 

We  thuB  find  v in  terms  of  n given  by  the  determinant:- 
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The  zero  order  approximation  gives  Aoo  c 
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6*  Inversion  of  Fourier  - Laplace  Transfonn. 


Bethe's  equation  16,  is,  in  our  notation 


J<  4i  V (z.u.n)  e 
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-:y  z 


(21) 


where  we  retain  the  variable  u instead  of  introducing  u*. 

We  may  now  proceed  to  solve  the  inhomogeneous  equation  (14)  at  the 
same  time  carrying  out  the  Bburier  inversion  of  $ as  in  Bethe’B  para  3. 
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and  A.  is  given  by 


dA(nv') 

rfl > 1 


To  obtain  the  Mellin  inversion  we  may  argue  as  follows: 

The  contribution  to  the  integral  in  (22)  will  arise  from  the 
neighbourhood  of  same  value  u0,  say,  of  u and  u0  will  be  a function 
of  n andz.  We  may  therefore  consider  n as  a function  of  uD  and  z 
so  that  using  (22)  we  have 


m -v(uo)z 

, , , , -Un(Uo)  e 

J du  % (zu)  e = 

0 go*(Uc)  ^ (uo) 


...  (25) 


Now  assume 
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where  (z,u),  go(u),  Mu)  are  all  slowly  varying  functions  of  u. 
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Then  we  may  write  (25)  in  the  fora 
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Hence 
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whence 

and  we  may  use  this  as  a 
general  relation  to  define  n for  a given  z and  u. 
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Now  expand  the  exponent  of  e in  (27)  in  powers  of  (u-uo) 

Thus  uhlu)  - n(uo}]  - z[v(u)  - v(uo)] 

(u-uo)9 
* 

2 

Htnce 


dn 

du 


d n 
+•  u — « 
du9 


da  v 
du9 


1 

X(z.uo) 


00 

”/■  du 
o 


(u-u0)9 

e 


dn  u 
du  + 8 


dan  _ z dV 
du4  2 du9 


✓ n 


d9n 

dav  . 

_ dn 

— 

“2  ) 

du9 

du9 

du 

F 

..  (27) 


..  (28) 


..  (29) 


..  (SO) 


. (31) 


7 
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But  we  have  from  (29)  u 
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so  that 
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Now  from  (29) 


d if 
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and  therefore 


Finally  we  have 


%(z,u) 


un  — zv 


Jin'*  A 


...  (34) 


Where  n is  given  by 
known  function  of  n 


u dv 
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for  given  u and  z since  i'  is  a 


Or 


4Hj(n+i) 


-4(n+i)  1 


A is  given  by  (23)  and  (24) 

Hence  the  result  is  not  altered  in  may  essential  way  by  taking 
into  account  capture  and  admixed  medium.  We  merely  get  a more  compli- 
cated expression  for  the  gs  and  hence  for  the  functions  v(n)  and  M"). 
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6.  Ihscussion 


'Hit;  chief  advantage  of  the  method  lies  in  its  ability  to  give 

T\ 

results  when  v and  are  not  small,  i.e.  at  distances  greater  than 
those  covered  by  age  theory  or  the  spherical  harmonic  method.^  However, 
it  is  interesting  to  apply  it  to  the  case  of  snail  v and  ff!  since 
the  other  methods  then  appear  as  successive  approximations.  ^ 


1.  Zero  order  approximation. 


This  is  given  by  equation  (00)  as 
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If  we  now  assume  that  v and 
lowest  powers  of  this  equation  we  have 
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are  small  for  all  j 


and  equate 
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for  large  Hj 


...  (35) 


We  have  2 cj  = c,  say  where  c is  the  number  of  scattered  neutrons 
J 

per  collision,  and  wc  may  write. 
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where  IT  is  the  average  inverse  mass  taken  over  the  scattering 
collisions. 

Thus 

v2  “'ll 

l-*--c-2cnfc  ...  (37) 
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We  must  therefore  assume  (l~c)  small  for  this  approximation 
i.e,  small  capture. 


From  (37) 
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Hence  the  approximation  is  valid  for  the  case  where  is  small 

i.e.  age  theory.  u 


UhT2V 

Now  consider  the  term  e in  equation  (34) 
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go  ~ c - 2 c If*  n which  we  may  take  as  unity  since  both  n H~  and 
(1  - c)  are  assumed  small. 

Thus  we  have 
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We  may  compare  this  expression  with  the  result  given  by  age 
velocity  theory  by  introducing  the  symbolic  age  t and  the  mean 
logarithmic  energy  loss  £ . Under  our  assumptions  with  respect  to  th<. 
cross  sections  we  find  that 


T a =7  + o(l)  • (45' 

Gc  bT 


£ - 8c  bf1  * O(jpr) 
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and  we  may  therefore  rewrite  equation  (44)  as 
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The  factor  exp 


is  the  usual  capture  correction 


- - (1  - c) 

obtained  by  standard  methods  (See  Marshak  op.  cit  pg.  821)  while  the 
remaining  factors  give  the  usual  age-velocity  solution  without  capture. 

Thus  the  age-velocity  solution  is  equivalent  to  our  equation  (44) 
when  we  take  account  of  the  dependence  of  t and  £ on  capture  and 
approximate  by  neglecting  terms  of  order/ in  their  expansions. 
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2.1st  order  approximation. 


Hie  first  order  approximation  to  equation  (1?)  is 


f(w) 


_1 

go 


If(v)  - 1] 


and  it  may  be  shown  that  this  equation  gives  n in  terms  of  v correctly 
up  to  order  ve  , If  we  expand  the  equation  up  to  terms  in  v and|  — j 

J 

we  get  the  following  expression  for  the  exponential  factor  in  (34) 
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where  p = (1  - c)  and  is  of  the  same  order  of  magnitude  as  v' 
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This  may  be  written,  neglecting  the  less  important  factors, 
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Harwell  Berks. 
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